General Properties of Two-dimensional Conformal 
Transformation in Electrostatics 



Yong Zeng/ Jinjie Liu,^ and Douglas H. Werner^ 
^Department of Electrical Engineering, 
Pennsylvania State University, University Park, PA 16802 
"^Department of Mathematical Sciences, 
Delaware State University, Dover, DE 19901 

Abstract 

Electrostatic properties of two-dimensional nanosystems can be described by their geometry 
resonances. In this paper we prove that these modes as well as the corresponding eigenvalues 
are invariant under any conformal transformation. This invariance further leads to a new way 
to studying the transformed structures. A special kind of three-dimensional transformations are 
further investigated. 



One specular property of Maxwell's equations is that their form will be invariant under ar- 
bitrary coordinate transformation if the field quantities as well as the material properties are 
transformed accordingly [l-3|, which leads to a powerful designing tool called transforma- 



tion optics or transformation electromagnetism in 2006 



2|]. Since then, a wealth of novel 



and unique devices have been theoretically suggested and/or experimentally fabricated, see 



recent reviews J^, [SJ and the references cited. An incomplete list includes e^ 
cloaks Q, y, M 



ectromagnetic 



ill 



event cloaks ^|, optical black holes 1(| and field splitters 
Quasi-static approximations can be applied when electromagnetic wavelength is far longer 
than the characteristic size of the material. The interactions between light and medium hence 
can be approximately described by Laplace's equation. Similar to the full-wave Maxwell 
equations, two-dimensional Laplacian equation is invariant under conformal transforma- 
tions 12| . Furthermore, the electrostatic transformation does not alter the properties of the 



constituent material, in sharp contrast to the full-wave coordinate transformations which 
generally transform normal medium to exotic and strange one. Lately, this unique Pfop- 
erty of Laplacian equation is used to study two-dimensional plasmonic particles ISHlSj. 
More specifically, two-dimensional complicated structures, by transforming to simple one- 
dimensional geometries such as finite-thickness metallic slab, can be studied analytically 



13 



18| . It is further found that these plasmonic nanoparticles, include crescent and touch- 



ing or non-touching cylinder dimers, can be broadband and significantly enhance the local 
electric fields. 

In this paper, we will study electrostatic conformal transformation in a new way, by taking 
advantage of the invariance of geometry resonances and their eigenvalues. We will show that 
transforming geometry can be equivalently explained as modifying the amplitudes of these 
geometry resonances by transforming the excitation source. Furthermore, a few energy 
quantities such as electrostatic energy possessed by the particle are found to be conserved 
under conformal transformation. Decreasing the area of the particle will then definitely 
increasing the degree of local field enhancement. 



We briefly recall the spectral Bergman-Milton theory 19|, l20|]. It is stated that the 
electrostatic behavior of a two-constituent composite can be described by a set of eigenmodes 
(fn of the following generalized eigenproblem 



V ■ [^(r)Vv?„,(r)] = s„VVn(r), 



(1) 



2 



where Sn representing the corresponding eigenvalues, and the function ^(r) characterizing 
the geometry of the composite, equal to 1 inside one constituent, with a permittivity of ei(w), 
and inside the other medium with a permittivity of e2(w) (we assume £2 = 1 in the following 
to simplify the discussion). Since this equation depends exclusively on the geometry, but 
not on the material composition, the resultant eigenmodes are therefore named as geometry 
resonances 19(]. Moreover, by defining a scalar product as 

(01^) = y V0* ■ dv, (2) 

we will find that Sn must be real and limited as < s„ < 1, and the normalized eigen- 
modes (fin form a complete orthonormal set [l9|. Note that this theory has been success fully 



employed to predict spaser 2l| and study ultrafast optical excitation in nanosystems 22|.|23|. 

It is important to notice, for a two-dimensional system, the eigenequation ([1]) remains 
invariance under any conformal coordinate transformation. To prove this fact, we rewrite 
the left-hand side of equation ([1]) as 

diBi ■ (edjifnej) = gijdi{9dj(pn) = gdi{9di(pn), (3) 

where the Einstein summation convention is employed, and the metric tensor Qij = ej • 
Sj = g6ij because of the conformal transformation [2]. Similarly the right-hand side can be 
reformulated as 

Sndiei ■ (djipnej) = Sngijdi{djlfn) = Sngduipn- (4) 

Evidently the transformed equation is identical to the previous one if we interpret the new 
equation as being in a right-handed Cartesian system and keep 6 and (pn unchanged. 

The invariance of the eigenmodes ipn and their corresponding eigenvalue s„ immediately 
suggest a way to studying the electrostatic response of the transformed structure. Note that 
we can expand its total potential (pt in terms of these eigenmodes, 

where k = u/c, s{uj) = 1/(1 — ci) being a spectral parameter and po being the external 



potential. The radiation damping has been included by adding the factor ik^/8 23, l25 |. 
Consequently, once these eigenmodes of the original or seed geometry are known, we can 
obtain the potential pt of any transformed geometry by just calculating the expansion co- 
efficients {ipn\<Po)- In other words, transforming geometry is equivalent to transforming the 



external source (po, or more precisely, changing the expansion coefficients of each eigenmodes. 
Furthermore, the time-averaged power absorbed by the structure can be written as 

= ^lm(ei) ^.„|/3„(v.„|y.o)r = ^Im(ei)/e, (6) 

n 

where /g represents the the integration of the electric field intensity \E\'^ inside the particle. 
In a similar way, the total power, the extinction, taken out by the particle has an expression 
as 

^ex = ^5^s„|(^„|^o)Plm[(6i - (7) 

n 

When we transform the external field (fo accordingly, i.e., (po{x^y) = (fQ{u,v) with u and v 
being the new coordinates, all these energy quantities. Pa? h and P^x, are evidently invariant. 
On the other hand, if the external potential is fixed in the w coordinate with w = u + iv,we 
will find that these energy quantities are proportional to when the conformal mapping 
has a form as a x w{z) with z = x + iy and a is real. In other words, bigger the particle 
area, stronger the absorption and extinction. 



It should 
the particle 



De mentioned that the geometry resonances actually are the surface modes of 



26 



27(1, and the resonance condition for the nth mode is strictly s{u) = Sn, 
or equivalently ei = 1 — when we do not include the radiation loss. Notice that the 

corresponding permittivity of the particle ei should be real and negative since < s„ < 1 



26| . The most striking property of the surface mode is that the resultant total electrostatic 



energy of the whole system, including the free space and the particle, is exactly zero 



28|. 



Furthermore, the complex frequency of the nth surface plasmonic resonance (SPR) is given 

by yiaii 

S{ujn - iln) = Sn, (8) 



with Un being real resonant frequency and 7„ being the relaxation rate [29|. Since the 
eigenvalues s„ are conserved under any conformal transformation, the transformed structure 
hence have same SPRs as the original structure at the identical resonant frequencies. For 
instance, since a metallic cylinder can be transformed from a metal-dielectric interface by 
w = e^, its SPR hence can be determined by nonretarded surface-plasmon condition, ei = 



■1, of the metal-dielectric interface 



3Q|. 



As a proof of principle, we consider an one- dimensional dielectric slab with finite thickness 
(Figure (la)) and its derivative systems obtained through different conformal transforma- 



^2 


^1 


(a), 


Q 
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FIG. 1: An one-dimensional finite-thickness slab in the xy coordinates (a) is transformed to an 
annulus (b) with a conformal transformation of w = e^, a crescent (c) and two kissing cylinders 
(d) with a conformal transformation of u; = l/z, in the new uv coordinates. Here w = u + iv and 
z = X + iy. 



tions. The eigenvalues of the slab are given by 

1 ih e^l^l('^2~'^i) 



Sk,± ^ , 

where k being real and nonzero, and the corresponding eigenmodes are 



-iky 



^g-2|fc|di -p g-|fc|(a!i+d2)]g|fc|x^ X <di 
^g-|/c|(di+d2)g|fc|x ^ g-|fc|x^ di<X <d2 

X > d2 



'1 ip ^\k\{d2-di)\^~\k\x 



(9) 



(10) 



where the normalize coefficient _|_ 



87r|fc|e"l^l'^i(e~l''l'^i =F e'^^^'^^). Evidently, in terms of 
the symmetry of electric field, these eigenmodes can be cataloged into two groups. The 
subscript +(— ) corresponds to even (odd) mode respectively, while their eigenvalues are 
bigger (smaller) than 0.5. Since the electric fields of the odd modes penetrate the metal 
weakly, they can propagate longer along the metallic surface than the even modes 30|. 
Furthermore, the resonant condition, approximately given by Re[s(a;„)] = s„ when the 



relaxation rate is weak, suggests that only the odd modes will be excited resonantly when 
the real part of the dielectric permittivity, Re(ei), is smaller than —1. Moreover, these 
eigenmodes and eigenvalues above can be directly applied to the complementary structure 
of the finite slab, i.e. a free-space gap sandwiched by two semi-infinite metallic spaces. 
Since these two structures have identical geometry, they therefore share same eigenvalues 
and eigenmodes. To obtain the induced potential of the complementary structure, we only 
need to interchange ei with €2 in the definition of s{oo). As a direct consequence, when 
Re(ei) is smaller than —1, only the even modes will be resonantly excited in the semi- infinite 



structure 



18|- 



We now calculate the resultant expansion coefficients of a structure transformed from the 
finite slab through a conformal mapping oiw = w{z). It is further assumed that the external 
potential ip^iu^v) = puU + PvV, which corresponding a uniform electric field —pu^u — Vv^v 
The corresponding expansion coefficients can be written as 

2k 

(V5fc,±lv5o) = [Q2pFk - giP*F^k] , (11) 

where 

f°° dw , 1 

dx dy—e''''', P=-{Pu-iPv)- (12) 

In addition, Qi = 1 and Q2 = ipe"'^'^"'^"'''^^^ for positive k, and they should exchange when 
A; < 0. 

To validate our approach, we first transforms the slab to two coaxial cylinders (Figure 
(lb)) with w = e^, and the radius of the two cylinders are ri = e'^^ and radius r2 = e'^'^ 
respectively. The expansion coefficients are found to be 

(^.,±1^0) = T^Pue'^^'^mi - 1), (13) 

here p^ is assumed to be zero because of the structural symmetry. Evidently, only eigenmodes 
with k = ±1 are excited. We further use these coefficients to calculate the induced potential 

f 

(rj - r|)(l - 62)^ r < ri 

r\{\ - 62)^ + r|(l - e\) - 2(1 - e2)r?ri/r^ n < r < r2 (14) 

rl{rl — rl){l — el)/r'^, r > r2 



with a = r|(l + €2)^ — rj(l — €2)"^- Our results are exactly identical to the one obtained 



by expanding the potential with Bessel functions 



26| . Furthermore, by letting ri — t- or 



ro — >■ oo, the results above can be used to describe a dielectric cylinder of a cylinder void 
26|. 

Since the eigenvalues of the finite slab cover a wide region [0, 1] (note that the eigenvalues 
of the eigenproblem ([T]) are limited between and 1), any negative Re(ei) will excite a surface 
mode as long as its expansion coefficient is not zero. Consequently the transformed system 
will present broadband response in principle. One example is the crescent studied in Ref. 
[isl . which can be obtained by using w = 1/z when rfi > (Figure (Ic)). The expansion 
coefficients are found to be 

Q/^ _|_ 

{^k,±\<fo) = \Vu + isgn(/c)p^] . (15) 

Evidently, each eigenmode is excited by the external potential. Furthermore, |(v5fc,±lv^o)| 
depends on the amplitude of the external electric field exclusively. The energy quantities 



such as the absorption hence do not depend on the direction of the incident field [l3|, a 



property which is not so obviously. Another example is the kissing cylinders suggested in 
Ref. [18|, which is obtain by transforming two semi-infinite slabs, with di < and d2 > 0, 
with w = 1/z (Figure (Id)). The corresponding expansion coefficients are calculated as 

(16) 



Q\k\{d2+d\) _j_ ^\k\{d2+di) -j- 

Pn-mi-jx^ + ^sgn(A;)p, 



' ^\k\{d2-di) I ' ""^'"^"'^'f'''g|fc|(d2-dl) ip 1 

Again all the eigenmodes can be excited by the external source. Note that when di = — ^2 
the cylinders have same radius. The kissing cylinders then possess both u and v mirror 
symmetry. Consequently, the even mode ipk,+ or the odd mode ipk- can be excited by u or 
V- polarized electric field alone HQ. ■ 

As suggested by the equation ([H]), the integration of the electric field intensity Jg inside 
the particle is invariant under conformal transformation. On the other hand, the area A 
of the particle is generally changed significantly. As a direct result, the degree of the local 
field localization, approximately measured by le/A, will be modified by the transformation. 
Taking the kissing cylinders or the crescent above as an example. Before transformation, the 
area of the one-dimensional slab is infinite. The transformed structure, on the other hand, 
has a finite area. Consequently, the kissing cylinders or the crescent strongly enhances the 



local field [li]l8\. 



A few words regarding three-dimensional transformation optics techniques in electro- 



static, which has been developed in Ref. 



3l| lately. Under a general transformation, the 



three-dimensional eigen-equation, equation ([T]), will be no longer invariant and the resultant 
expression usually is very complicated. One exception is that the Jacobian matrix J of the 
transformation satisfies J = f^^{r')X, with X being the unit 3x3 matrix and f{r') being 
an arbitrary function of r' (the transformation used in Ref. 3l| possesses this property). 
Under this condition, the new eigenmodes 0„ connect to the previous eigenmodes ipn of the 
original system as 

V</.„ = y^V^n. (17) 
The corresponding eigenvalues t„ are also different with the previous Sn 

when the integrations are performed over the previous system. Once the new eigenmodes and 
eigenvalues are obtained, equation (5-7) can be directly employed without any modification. 
Note that for a two-dimensional conformal transformation equation (17,18) still work, by 
simply replacing the function / with a constant number. Again, we find that the eigenmodes 
and eigenvalues are invariant under any two-dimensional conformal transformation. 

In conclusion, we proved that the electrostatic eigenmodes and eigenvalues of two- 
dimensional nanosystems are invariant under any conformal transformation. Based on this 
property, we suggested an new approach to studying the electrostatic responses of the trans- 
formed structures. Namely, transforming a geometry is equivalent to transforming the exter- 
nal potential or the expansion coefficients of the invariant eigenmodes. Furthermore, except 
the potential, energy quantities such as absorption, extinction and electrostatic energy are 
also found to be conserved. Moreover, the significant conditions to designing broadband 
nanosystems are the eigenvalues of the system should cover a wide portion of [0, 1] as well 
as the external source should excite nearly all the eigenmodes. 

This work was supported in part by the Penn State MRSEC under NSF grant no. DMR 
0213623. 
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